Probing Hard Color- Singlet Exchange in pp Collisions 
at v/i = 630 GeV and 1800 GeV 



B. Abbott, 40 M. Abolins, 37 V. Abramov, 15 B.S. Acharya, 8 I. Adam, 39 D.L. Adams, 48 
M. Adams, 24 S. Ami, 23 H. Aihara, 17 G.A. Alves, 2 N. Amos, 36 E.W. Anderson, 30 R. Astur, 42 
M.M. Baarmand, 42 V.V. Babintsev, 15 L. Babukhadia, 16 A. Baden, 33 B. Baldin, 23 
S. Banerjee, 8 J. Bantly, 45 E. Barberis, 17 P. Baringer, 31 J.F. Bartlett, 23 A. Belyaev, 14 
S.B. Beri, 6 I. Bertram, 26 V.A. Bezzubov, 15 P.C. Bhat, 23 V. Bhatnagar, 6 
M. Bhattacharjee, 42 N. Biswas, 28 G. Blazey, 25 S. Blessing, 21 P. Bloom, 18 A. Boehnlein, 23 
N.I. Bojko, 15 F. Borcherding, 23 C. Boswell, 20 A. Brandt, 23 R. Breedon, 18 R. Brock, 37 
A. Bross, 23 D. Buchholz, 26 V.S. Burtovoi, 15 J.M. Butler, 34 W. Carvalho, 2 D. Casey, 37 
Z. Casilum, 42 H. Castilla-Valdez, 11 D. Chakraborty, 42 S.-M. Chang, 35 S.V. Chekulaev, 15 
W. Chen, 42 S. Choi, 10 S. Chopra, 36 B.C. Choudhary, 20 J.H. Christenson, 23 M. Chung, 24 
D. Claes, 38 A.R. Clark, 17 W.G. Cobau, 33 J. Cochran, 20 L. Coney, 28 W.E. Cooper, 23 
C. Cretsinger, 41 D. Cullen-Vidal, 45 M.A.C. Cummings, 25 D. Cutts, 45 O.I. Dahl, 17 
K. Davis, 16 K. De, 46 K. Del Signore, 36 M. Demarteau, 23 D. Denisov, 23 S.P. Denisov, 15 
H.T. Diehl, 23 M. Diesburg, 23 G. Di Loreto, 37 P. Draper, 46 Y. Ducros, 5 L.V. Dudko, 14 
\D '. S.R. Dugad, 8 A. Dyshkant, 15 D. Edmunds, 37 J. Ellison, 20 V.D. Elvira, 42 R. Engelmann, 42 

S. Eno, 33 G. Eppley, 48 P. Ermolov, 14 O.V. Eroshin, 15 V.N. Evdokimov, 15 T. Fahland, 19 
M.K. Fatyga, 41 S. Feher, 23 D. Fein, 16 T. Ferbel, 41 G. Finocchiaro, 42 H.E. Fisk, 23 
Y. Fisyak, 43 E. Flattum, 23 G.E. Forden, 16 M. Fortner, 25 K.C. Frame, 37 S. Fuess, 23 
E. Gallas, 46 A.N. Galyaev, 15 P. Gartung, 20 V. Gavrilov, 13 T.L. Geld, 37 R.J. Genik II, 37 
K. Genser, 23 C.E. Gerber, 23 Y. Gershtein, 13 B. Gibbard, 43 B. Gobbi, 26 B. Gomez, 4 
G. Gomez, 33 P.I. Goncharov, 15 J.L. Gonzalez Soli's, 11 H. Gordon, 43 L.T. Goss, 47 
K Gounder, 20 A. Goussiou, 42 N. Graf, 43 PD. Grannis, 42 D.R. Green, 23 H. Greenlee, 23 
instein, 1 P. Grudberg, 17 S. Griinendahl, 23 G. Guglielmo, 44 J.A. Guida, 16 J.M. Guida, 45 
lupta, 8 S.N. Gurzhiev, 15 G. Gutierrez, 23 P. Gutierrez, 44 N.J. Hadley, 33 H. Haggerty, 23 



J.M 



S. Grinstein, 1 P. Grudberg, 17 S. Griinendahl, 23 G. Guglielmo, 44 J.A. Guida, 16 J.M. Guida 

A. Gupta, 8 S.N. Gurzhiev, 15 G. Gutierrez, 23 P. Gutierrez, 44 N.J. Hadley, 33 H. Haggerty, 
S. Hagopian, 21 V. Hagopian, 21 K.S. Hahn, 41 R.E. Hall, 19 P. Hanlet, 35 S. Hansen, 23 

J.M. Hauptman, 30 D. Hedin, 25 A. P. Heinson, 20 U. Heintz, 23 R. Hernandez-Montoya, 11 
T. Heuring, 21 R. Hirosky, 24 J.D. Hobbs, 42 B. Hoeneisen, 4 '* J.S. Hoftun, 45 F. Hsieh, 36 
Ting Hu, 42 Tong Hu, 27 A.S. Ito, 23 E. James, 16 J. Jaques, 28 S.A. Jerger, 37 R. Jesik, 27 
T. Joffe-Minor, 26 K. Johns, 16 M. Johnson, 23 A. Jonckheere, 23 M. Jones, 22 H. Jostlein, 23 
S.Y. Jun, 26 C.K. Jung, 42 S. Kahn, 43 G. Kalbfleisch, 44 D. Karmanov, 14 D. Karmgard, 21 
R. Kehoe, 28 M.L. Kelly, 28 S.K. Kim, 10 B. Klima, 23 C. Klopfenstein, 18 W. Ko, 18 
J.M. Kohli, 6 D. Koltick, 29 A.V. Kostritskiy, 15 J. Kotcher, 43 A.V. Kotwal, 39 A.V. Kozelov, 15 
E.A. Kozlovsky, 15 J. Krane, 38 M.R. Krishnaswamy, 8 S. Krzywdzinski, 23 S. Kuleshov, 13 
Y. Kulik, 42 S. Kunori, 33 F. Landry, 37 G. Landsberg, 45 B. Lauer, 30 A. Leflat, 14 J. Li, 46 

Q.Z. Li-Demarteau, 23 J.G.R. Lima, 3 D. Lincoln, 23 S.L. Linn, 21 J. Linnemann, 37 
R. Lipton, 23 F. Lobkowicz, 41 S.C. Loken, 17 A. Lucotte, 42 L. Lueking, 23 A.L. Lyon, 33 
A.K.A. Maciel, 2 R.J. Madaras, 17 R. Madden, 21 L. Magana-Mendoza, 11 V. Manankov, 14 
S. Mani, 18 H.S. Mao, 23 '* R. Markeloff, 25 T. Marshall, 27 M.I. Martin, 23 K.M. Mauritz, 30 

B. May, 26 A.A. Mayorov, 15 R. McCarthy, 42 J. McDonald, 21 T. McKibben, 24 J. McKinley, 37 
T. McMahon, 44 H.L. Melanson, 23 M. Merkin, 14 K.W. Merritt, 23 C. Miao, 45 H. Miettinen, 48 
A. Mincer, 40 C.S. Mishra, 23 N. Mokhov, 23 N.K. Mondal, 8 H.E. Montgomery, 23 P. Mooney, 4 



1 



M. Mostafa, 1 H. da Motta, 2 C. Murphy, 24 F. Nang, 16 M. Narain, 23 V.S. Narasimham, 8 

A. Narayanan, 16 H.A. Neal, 36 J.P. Negret, 4 P. Nemethy, 40 D. Norman, 47 L. Oesch, 36 
V. Oguri, 3 E. Oliveira, 2 E. Oltman, 17 N. Oshima, 23 D. Owen, 37 P. Padley, 48 A. Para, 23 

Y.M. Park, 9 R. Partridge, 45 N. Parua, 8 M. Paterno, 41 B. Pawlik, 12 J. Perkins, 46 
M. Peters, 22 R. Piegaia, 1 H. Piekarz, 21 Y. Pischalnikov, 29 B.G. Pope, 37 H.B. Prosper, 21 
S. Protopopescu, 43 J. Qian, 36 P.Z. Quintas, 23 R. Raja, 23 S. Rajagopalan, 43 O. Ramirez, 24 

S. Reucroft, 35 M. Rijssenbeek, 42 T. Rockwell, 37 M. Roco, 23 P. Rubinov, 26 R. Ruchti, 28 
J. Rutherfoord, 16 A. Sanchez-Hernandez, 11 A. Santoro, 2 L. Sawyer, 32 R.D. Schamberger, 4 
H. Schellman, 26 J. Sculli, 40 E. Shabalina, 14 C. Shaffer, 21 H.C. Shankar, 8 R.K. Shivpuri, 7 
D. Shpakov, 42 M. Shupe, 16 H. Singh, 20 J.B. Singh, 6 V. Sirotenko, 25 E. Smith, 44 
R.P Smith, 23 R. Snihur, 26 G.R. Snow, 38 J. Snow, 44 S. Snyder, 43 J. Solomon, 24 
M. Sosebee, 46 N. Sotnikova, 14 M. Souza, 2 G. Steinbriick, 44 R.W. Stephens, 46 
M.L. Stevenson, 17 D. Stewart, 36 F. Stichelbaut, 42 D. Stoker, 19 V. Stolin, 13 
D.A. Stoyanova, 15 M. Strauss, 44 K. Streets, 40 M. Strovink, 17 A. Sznajder, 2 P. Tamburello; 
J. Tarazi, 19 M. Tartaglia, 23 T.L.T. Thomas, 26 J. Thompson, 33 T.G. Trippe, 17 P.M. Tuts, 3 
V. Vaniev, 15 N. Varelas, 24 E.W. Varnes, 17 D. Vititoe, 16 A. A. Volkov, 15 A.P Vorobiev, 15 
H.D. Wahl, 21 G. Wang, 21 J. Warchol, 28 G. Watts, 45 M. Wayne, 28 H. Weerts, 37 A. White, 4 

J.T. White, 47 J.A. Wightman, 30 S. Willis, 25 S.J. Wimpenny, 20 J.V.D. Wirjawan, 47 
J. Womersley, 23 E. Won, 41 D.R. Wood, 35 Z. Wu, 23 ^ R. Yamada, 23 P. Yamin, 43 T. Yasuda, : 
P. Yepes, 48 K. Yip, 23 C. Yoshikawa, 22 S. Youssef, 21 J. Yu, 23 Y. Yu, 10 B. Zhang, 23 <t 
Y. Zhou, 23 -* Z. Zhou, 30 Z.H. Zhu, 41 M. Zielinski, 41 D. Zieminska, 27 A. Zieminski, 27 

E.G. Zverev, 14 and A. Zylberstejn 5 



(D0 Collaboration) 



1 Universidad de Buenos Aires, Buenos Aires, Argentina 
2 LAFEX, Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, Brazil 
3 Universidade do Estado do Rio de Janeiro, Rio de Janeiro, Brazil 

4 Universidad de los Andes, Bogota, Colombia 
5 DAPNI A/Service de Physique des Particules, CEA, Saclay, France 
e Panjab University, Chandigarh, India 
7 Delhi University, Delhi, India 
8 Tata Institute of Fundamental Research, Mumbai, India 
9 Kyungsung University, Pusan, Korea 
10 Seoul National University, Seoul, Korea 
n CINVESTAV, Mexico City, Mexico 
12 Institute of Nuclear Physics, Krakow, Poland 
13 Institute for Theoretical and Experimental Physics, Moscow, Russia 
14l Moscow State University, Moscow, Russia 
15 Institute for High Energy Physics, Protvino, Russia 
16 University of Arizona, Tucson, Arizona 85721 
17 Lawrence Berkeley National Laboratory and University of California, Berkeley, California 94720 

18 University of California, Davis, California 95616 

19 University of California, Irvine, California 92697 



2 



20 University of California, Riverside, California 92521 

21 Florida State University, Tallahassee, Florida 32306 

22 University of Hawaii, Honolulu, Hawaii 96822 
23 Fermi National Accelerator Laboratory, Batavia, Illinois 60510 
24 University of Illinois at Chicago, Chicago, Illinois 60607 
25 Northern Illinois University, DeKalb, Illinois 60115 
26 Northwestern University, Evanston, Illinois 60208 
27 Indiana University, Bloomington, Indiana 47405 
28 University of Notre Dame, Notre Dame, Indiana 46556 
29 Purdue University, West Lafayette, Indiana 47907 
30 Iowa State University, Ames, Iowa 50011 
31 University of Kansas, Lawrence, Kansas 66045 
32 Louisiana Tech University, Ruston, Louisiana 71272 
33 University of Maryland, College Park, Maryland 20742 

34 Boston University, Boston, Massachusetts 02215 
35 Northeastern University, Boston, Massachusetts 02115 
36 University of Michigan, Ann Arbor, Michigan 48109 
37 Michigan State University, East Lansing, Michigan 48824 

38 University of Nebraska, Lincoln, Nebraska 68588 

39 Columbia University, New York, New York 10027 
A0 New York University, New York, New York 10003 

41 University of Rochester, Rochester, New York 14627 
42 State University of New York, Stony Brook, New York 11794 
43 Brookhaven National Laboratory, Upton, New York 11973 

44 University of Oklahoma, Norman, Oklahoma 73019 

45 Brown University, Providence, Rhode Island 02912 

46 University of Texas, Arlington, Texas 76019 
47 Texas A&M University, College Station, Texas 77843 
48 Rice University, Houston, Texas 77005 



Abstract 



We present results on dijet production via hard color-singlet exchange in 
proton-antiproton collisions at ^/s = 630 GeV and 1800 GeV using the D0 
detector. The fraction of dijet events produced via color-singlet exchange is 
measured as a function of jet transverse energy, separation in pseudorapidity 
between the two highest transverse energy jets, and proton-antiproton center- 
of-mass energy. The results are consistent with a color-singlet fraction that 
increases with an increasing fraction of quark-initiated processes and incon- 
sistent with two-gluon models for the hard color-singlet. 
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The exchange of a quark or gluon between interacting partons in hadronic collisions 
typically results in final-state particle production over several units of rapidity In contrast, 
the exchange of a color-singlet is expected to yield a rapidity gap, defined as the absence of 
particles in a region of rapidity (or pseudorapidity, rj = — lntan(#/2), where 9 is the polar 
angle). Rapidity gaps are most commonly observed in low momentum-transfer diffractive 
and elastic scattering processes, which are attributed to the exchange of a color-singlet called 
the pomeron. The observation of rapidity gaps between jets at both the Fermilab pp collider 
(Tevatron) |^@| arid the DESY ep Collider (HERA) || implies the exchange of a hard 
color-singlet. The measured fraction of dijet events arising from color-singlet exchange is 
roughly 1% in proton-antiproton collisions []1|-§J and 10% in positron-proton collisions ||. 
These rates are too large to be explained by electroweak boson exchange and indicate a 
strong-interaction process [0,0. 

The final-state particle multiplicity distribution between jets is used to distinguish be- 
tween color exchange and color-singlet exchange. By identifying an excess of low multi- 
plicity events using the D0 calorimeter, tracking system, or both, we measure the observ- 
able fraction of color-singlet exchange in dijet events ("color-singlet fraction"). This Letter 
presents new measurements by the D0 Collaboration of the color-singlet fraction as a func- 
tion of jet transverse energy (E T ), pseudorapidity separation between the two highest E T 
jets (A77 = \rji — 772 1), and proton-antiproton center-of-mass energy (y/s). We then compare 
these measurements to hard color-singlet models. 

The D0 detector and trigger system are described in Ref. ||. The data samples were 
accumulated during the 1994-1996 Tevatron run primarily at a proton-antiproton center-of- 
mass energy of 1800 GeV, with a short run at 630 GeV during this period. Jets were found 
in the liquid-argon calorimeter (full coverage to \rj\ <4.1) using a cone algorithm with radius 
TZ— 0.7 in rj — (ft space, where is the azimuthal angle. Samples of "opposite-side" (OS) jet 
events (771 ■ r] 2 < 0) and "same-side" (SS) jet events (7/1 • 772 > 0) were collected using triggers 
similar to those described in Ref. 0. The triggers required both jets to have E T above a 
threshold which varied from 12 to 25 GeV depending on the trigger. All triggers required at 
least two jets with \r]\ > 1.6, and the opposite-side triggers required a dijet pseudorapidity 
separation of Ar/>4.0 at y/s = 1800 GeV and Ar/>3.2 at 630 GeV. 

In the offline analysis, the vertex was required to have a longitudinal position within 
50 cm of the center of the detector and the two highest transverse energy (leading) jets were 
required to have \r)\ > 1.9. The opposite-side samples were required to have A7/>4.0 between 
the leading two jets. Since multiple interactions (more than one proton-antiproton inter- 
action in the same bunch crossing) may obscure the rapidity gap signal of color-singlet ex- 
change, we retained only events with a single reconstructed vertex and additionally required 
the beam-beam hodoscope timing information to be consistent with a single interaction. 

Table | summarizes the data samples used in this analysis, including the cut on the trans- 
verse energy of the second leading jet (£^2), the integrated luminosity, and the approximate 
number of events for each trigger. At y/s = 1800 GeV, three exclusive opposite-side data 
samples corresponding to triggers with low, medium, and high Et thresholds ("low-Er", 
"med-£y\ and "high-£y') are used to determine the E T dependence of the color-singlet 
fraction. The A77 dependence of the color-singlet fraction is measured using the low-Er and 
high-E^ samples. The large statistics high-E'r sample gives the most accurate measure of 
the color-singlet fraction. 
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Samples with the same jet Et threshold of 12 GeV are used to compare the = 630 GeV 
and 1800 GeV color-singlet fractions. This lower E T2 cut provides adequate statistics for 
the 630 GeV sample. Figures |l](a) and (b) show the average jet Et and Ar/, respectively, 
of the two leading jets for the 630-OS and 1800-OS samples. As expected from kinematic 
considerations, the jets in the \fs = 1800 GeV sample reach a higher Et and Arj than the 
630 GeV sample. Figure |l](c) shows the average parton x (x) for the two samples, where we 
calculate x using a leading order approximation (x = (x + +x^)/2; x± = Yui=i ^Ti^ ±r>1 / \/s) ■ 

Table [I] also lists the same-side samples at 630 GeV and 1800 GeV. Same-side events are 
useful since they have similar characteristics to opposite-side events but have no observable 
color-singlet component ||. 



Name 


(GeV) 


E T 2 (GeV) 


Luminosity (nb 1 ) 


Events 


low-Er 


1800 


15 - 25 


86 


27000 


med-ET 




25 - 30 


1300 


21000 


high-Er 




> 30 


11000 


72000 


1800-OS 




> 12 


86 


48000 


1800-SS 




> 12 


18 


57000 


630-OS 


630 


> 12 


520 


6700 


630-SS 




> 12 


350 


6800 



TABLE I. Data samples showing the offline cut on the transverse energy of the second leading 
jet (£t2)> integrated luminosity, and approximate number of events after all offline cuts. All triggers 
are opposite-side jet triggers except for 1800-SS and 630-SS, which are same-side jet triggers. 




FIG. 1. Jet characteristics of the 630-OS (solid line) and 1800-OS (dashed line) data samples. 
The normalized distributions are shown for (a) the average Et, (b) A 77, and (c) x of the two leading 
jets. 

We determine the amount of color-singlet exchange in the opposite-side samples by uti- 
lizing the multiplicity distributions measured by the calorimeter and the tracking system 
in the pseudorapidity region \r)\ < 1. This region has little sensitivity to contamination 
from the leading jets, which have a minimum \rj\ of 1.9. We measure the multiplicity in 
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the electromagnetic calorimeter by counting the number of towers (At/ x A0 = 0.1 x 0.1) 
with Et>200 MeV (n C ai) ■ The electromagnetic section of the calorimeter has a low level of 
noise and the ability to detect both charged and neutral particles. We also count the num- 
ber of tracks in the central drift chamber (n tr ^), which is a non- magnetic tracking system 
that detects charged particles to low momenta and provides an independent measure of the 
multiplicity. 

Figure |^(a) displays the Ti ca i versus ra tr k distribution in the central region (|t/| < 1) between 
the jets for the high-i^ (y/s = 1800 GeV) sample. A large excess of events is observed at low 
multiplicity, consistent with the expectations of a color-singlet exchange signal and similar 
to that observed previously [0,0. Figure 0(b) shows the n ca \ distribution with a three- 
parameter negative binomial distribution (NBD) fit, which is used to determine the color- 
exchange background. NBD's have previously been shown to provide a good description of 
particle multiplicity distributions in proton-antiproton collisions ||. 




FIG. 2. Multiplicity in the region \rj\ < 1 between the two leading jets for the high-E'r sample: 
(a) two-dimensional multiplicity, n ca \ vs. nt r k! (b) n ca \ only with NBD fit. 

We measure the color-singlet fraction (fg) by dividing the number of events above the 
NBD fit for n ca i < 1 by the total number of events in the sample. This quantity is insensitive 
to noise: after the suppression of a few known noisy calorimeter towers, less than 1% of 
events taken during no-beam runs have more than one tower above threshold. In addition, 
the color-singlet fraction is not sensitive to jet finding and trigger efficiencies, which cancel 
in the ratio since the jets found in color-singlet events are indistinguishable from the jets in 
color-exchange events. 

To minimize sensitivity to residual contamination from multiple proton-antiproton inter- 
actions, we fit to the low multiplicity region of the data. The fits start where the signal to 
background ratio becomes negligible (n ca \ = 2) and end near the maximum of the multiplicity 
distribution (n ca i = 11 at 630 GeV and n ca i = 14 at 1800 GeV). The fits are then extrapolated 
to zero multiplicity to determine the background. This method was tested by generating 
and fitting ensembles of multiplicity distributions with and without a known signal compo- 
nent at low multiplicity. These studies confirm that the fitting method is unbiased and also 
provide a measure of the systematic error J7^Q . An additional uncertainty due to the choice 
of the starting multiplicity bin used for the fit is estimated by determining the background 
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using a starting bin of n ca i = 3; the resulting change in background extrapolation is applied 
as an additional systematic error. 

Corrections are applied to f$ to account for residual multiple interaction contamination 
in the samples as well as for single- vertex events erroneously tagged as multiple interaction 
events. The extra multiplicity from secondary interactions causes color-singlet events to be 
erroneously measured with a large multiplicity. By studying the shape of the multiplicity 
distribution as a function of luminosity, a correction is determined that results in roughly 
a 30% increase in f$. The vertex correction reduces fs by about 10% since color-singlet 
events always have one vertex, but about 10% of color exchange events with one interaction 
have two reconstructed vertices [[|. 

Applying the fitting method to the = 1800 GeV high-£ T data (E T2 > 30 GeV) 
shown in Fig. 0(b), we obtain f s = [0.94 ± 0.04(stat)±0.12(syst)]%. The systematic error 
is dominated by the uncertainty in the background subtraction (0.09%) but also includes 
uncertainties in the single interaction determination (0.07%), offline cuts (0.03%), and jet 
energy scale and reconstruction (0.02%). This measurement is the most precise to date and 
is consistent with previously published measurements at the Fermilab Tevatron with similar 
minimum jet E? !]-§]]• 

The value of fs does not depend on the exact method used to measure particle multi- 
plicity. We vary the calorimeter tower threshold used for determining n ca i between 150 MeV 
and 350 MeV (the nominal value is 200 MeV) and obtain a consistent value for the fraction. 
Redefining multiplicity using a "cluster" of neighboring calorimeter towers to account for 
possible calorimeter showering effects, or using the tracking multiplicity (n tr k) instead of n ca \ 
results in a measured color-singlet fraction within a few per cent of the nominal value. 

Figure |3|(a) shows the multiplicity distribution between the jets for the 1800-OS sample 
with Et2 > 12 GeV. Fitting the n ca \ distribution and measuring the color-singlet fraction, we 
obtain fs = [0.54 ± 0.06(stat)±0.16(syst)]%. This value is about 2a smaller than the value 
measured for E?2 >30 GeV, indicating a dependence of the color-singlet fraction on jet E?- 

Figure |(b) shows the multiplicity distribution between the jets for the 630-OS sample. 
For this sample we obtain fs = [1.85 ± 0.09(stat)±0.37(syst)]%. The uncertainties in the 
630-OS and 1800-OS measurements are dominated by the fit errors, which are not correlated 
between the two samples. We calculate the significance of the difference between the two 
values as 3.1<j in the limiting case that the uncertainties are completely uncorrelated. 




FIG. 3. Two-dimensional multiplicity (n ca i vs. n tr k) in the region \r/\ < 1 for the (a) 1800-OS 
and (b) 630-OS samples. 
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The ratio of the two measurements for E T2 > 12 GeV is Rf 8 ° 00 = / s (630)// s (1800) = 
3.4 ± 1.2. We verify that the lower average multiplicity for 630 GeV does not affect the 
measured ratio by raising the electromagnetic tower threshold of the 1800 GeV data so that 
the leading edge of the n cal distribution matches that of the 630 GeV sample. We also obtain 
a consistent value for the ratio when we measure it using a two-dimensional n c& \ versus n tr k 
multiplicity method (described below). 

Measuring the color-singlet fraction as a function of Et and Arj requires splitting the 
data into several bins. The NBD fitting uncertainty becomes unacceptably large for these 
smaller statistics sub-samples, so a new two-dimensional multiplicity method is employed to 
reduce the uncertainty. For each bin of Ej- or Ar/, we first determine the fraction of events 
with n ca i + 72 t rk < 1- We then obtain the "2D" color-singlet fraction (f^^) by subtracting 
the appropriate color-exchange background and correcting by the acceptance relative to the 
NBD method. 

The color-exchange background for the "2D" fraction is determined using the following 
procedure. The NBD fit to the calorimeter multiplicity distribution for each full sample is 
used to determine the total background (n. ca i<l) for that sample. This "ID" background is 
converted to a "2D" background by multiplying by the ratio N(n ca \ + n tr k < 1)/N{n ca \ < 1) 
obtained from the same-side multiplicity distribution in the same detector region. The 
same-side data are thus used only to account for correlations between n cai and n trk . The 
average background is roughly 25% for each sample. The dependence of the background on 
Et and Ar? is measured using a low-multiplicity control region (3<rr ca i + rr tr k<5) where the 
color-singlet contribution is negligible. The final background for each sample is determined 
by multiplying the average background with its Et and At? dependence. 

The background is then subtracted from each Et and Ar] bin of the uncorrected "2D" 
fraction. The background has little dependence on E T and Ar/, so the dependence of the 
color-singlet fraction on these variables before and after background subtraction is quite 
similar. 

Comparing the color-singlet fraction from the "2D" method to fs from the NBD fitting 
method for each sample, we determine that the "2D" method has a relative acceptance of 
80% for color-singlet events. This acceptance has been determined to have little dependence 
on Et and Arj, and its uncertainty (3% over the Et range and 6% over the Arj range) 
is included in subsequent shape comparisons to different color-singlet models. We obtain 
the final color-singlet fractions fg (Et) and fj? D \Ar]) by correcting for the color-singlet 
acceptance determined from the high-i?^ sample. 

Figure ^(a) shows the color-sing let fraction ff D) (E T ) at 1800 GeV obtained using the 
\ow-Et, med-.&r, and high-ET data samples. The error bars show the statistical uncertainty, 
while the bands show the relative normalization uncertainties between the different data 
samples. We bin by second leading jet Et (Et2) to mimic the way the dijet triggers select 
events. In contrast, binning by average dijet Et selects primarily back-to-back dijet events 
near the trigger Et threshold. Since color-singlet events are more likely to be back-to-back 
than color-exchange events, using an average Et method would bias the color-singlet fraction 
as a function of Et- 

Figures |](b) and f|(c) show the color-singlet fraction ff D \Arf) at 1800 GeV for the low- 
Et and high-.&r samples, respectively. The error bars show the statistical uncertainty, while 
the bands show the relative normalization uncertainties between the different data samples. 
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Figures |](a)-(c) show that the measured color-singlet fraction tends to increase with Et 
and At], implying that the fraction increases with parton x. Replotting Fig. £|(b)-(c) in 
Fig. f|(d) as a function of x, with x averaged over each At] bin, shows a rise in the measured 
color-singlet fraction with x. In principle it is possible to measure the color-singlet fraction 
directly as a function of parton x. However, bin-dependent acceptance effects related to the 
dijet trigger make it preferable to do the analysis in the observables Et and Ai]. 
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FIG. 4. The color-singlet fraction f^ D ^: (a) as a function of the second leading jet Et; as a 
function of Arj between the two leading jets for (b) the low- 1st sample and (c) the high-i^T sample; 
(d) as a function of x for each Arj bin in (b) and (c). The horizontal error bar on the first x 
bin shows the la spread in x values entering that (representative) bin. Statistical error bars and 
relative normalization uncertainties for each sample (hatched bands) are shown. 
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Measuring the color-singlet fraction as a function of Et, At], and yfs probes the nature 
of hard color-singlet exchange. If the color-singlet exchange process were identical to single- 
gluon exchange except for different coupling factors to quarks and gluons, the color-singlet 
fraction would depend only on parton distribution functions through the parton x. Thus, 
for a color-singlet that couples more strongly to gluons (quarks), the color-singlet fraction 
would fall (rise) as a function of increasing x, since the gluon distribution becomes relatively 
suppressed as x increases. This scenario would imply a decreasing (increasing) color-singlet 
fraction with increasing jet Et and Arj or decreasing yfs. The observed color-singlet fraction 
may deviate from this simplistic behavior, however, if the dynamics of color-singlet exchange 
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are distinct from the dynamics of single-gluon exchange. Currently QCD does not account 
for the existence of hard color-singlet exchange, but higher-order QCD processes have been 



proposed to explain this phenomenon fll0|-|T6| . 

The measured color-singlet fraction includes the probability that the color-singlet events 
are observable, and can be written as fs — S ■ Os ing i ct /a, where cr sing i et /a is the fraction of 
dijet events produced by color-singlet exchange, and S is the survival probability (proba- 
bility that the rapidity gap is not contaminated by particles from spectator interactions). 
This survival probability (5*^10%) cannot be measured independently from the color sin- 
glet, and is assumed not to depend on the x or type of the initiating partons in the hard 



scattering |fTT]JT7[1 . Although initial estimates of 5* had a weak (~ 1/lns) dependence on 



center-of-mass energy |l]J , a recently revised estimate gives a stronger dependence (18 



We compare our color-singlet results to Monte Carlo color-singlet models using HER- 



WIG ||19|| , which includes parton showering and hadronization of the final state partons. 
Monte Carlo samples were generated using CTEQ2M parton distribution functions |f2T[ . 
As in the data, jets were reconstructed using a cone algorithm with a radius of 1Z = 0.7. 
Detector effects were simulated by smearing the jet Et and rj according to measured jet 
resolutions. The two leading jets were required to have \r]\ > 1.9 and Arj > 4.0 as in the 
data samples. 

For all the models the color-singlet fraction as a function of Et and Ar] is obtained 
by dividing the binned color-singlet cross section by the corresponding HERWIG QCD cross 
section. A rapidity gap requirement was not imposed on any of the Monte Carlo samples 
since identifying rapidity gaps depends on the details of hadronization and the modeling of 
spectator interactions. The rapidity gap survival probability is assumed to be independent 
of Et and Ar], and its magnitude is absorbed in an overall normalization factor. 

We use a Bayesian probability method to fit the color-singlet fraction predicted by the 
Monte Carlo models independently to the J^ D \Et) data (Fig. £|(a)) and the ff D \Arj) 
data (Fig. |](b)-(c)). Although the measurements of the color-singlet fraction as a function 
of Et and Arj are statistically correlated, they contain complementary information. The 
normalization, which is not well-determined for any of the models, is the only free parameter 
in the fits to each model, except for the "free- factor" model (discussed below). 

The corrections, background, acceptance, Monte Carlo expectations, and uncertainties in 
these quantities are included in the Bayesian integral as prior probability distributions. Since 
there has been no previous comparison to color-singlet models, we assign an identical prior 
probability for each model and assume a flat prior probability (no preferred value) for each 
fit parameter. We obtain the posterior probabilities for the fit parameters by integrating the 
likelihood function over all prior probability distributions. For each model, we determine the 
fit parameter values that maximize the posterior probability for that model. The uncertainty 
in each fit parameter is obtained from its posterior probability distribution by determining 
the 68% confidence interval which has equally probable endpoints. If one of the boundaries 
is restricted by physical constraints then we use the smallest 68% confidence interval to 
determine the uncertainty. The maximum of the likelihood is taken as a measure of the 
relative goodness-of-fit of that model. 

The Bayesian method was chosen for fitting the data since the conventional \ 2 method 



introduces a normalization bias in the case of correlated uncertainties 22 . This bias is small 



when the model describes the data well, but for some of the models the fitted normalization 
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from the x 2 method was as much as a factor of two smaller than the corresponding Bayesian 
value. 

Figures [5]-|6| show fits to the measured color-singlet fraction. Even though the fraction 
versus Ar] and E? are fit separately, the fit parameters obtained from each variable do not 
differ appreciably for any given model. Since the Ar] distributions for the low--Er and high- 
Et samples are fit simultaneously, a model prediction that is above the data for one plot 
would be below the data for the other. 

The dashed line in Fig. [5] (BFKL 2^2) shows a parton-level color-singlet prediction of a 
two-gluon exchange process incorporating the BFKL dynamics of Mueller and Tang [13, 23]. 



The solid line is the BFKL color-singlet fraction at the jet level after hadronization. The 
BFKL models do not give an adequate description of the data, primarily due to their predic- 
tion of decreasing color-singlet fractions with Et- Assuming a survival probability of about 
10%, the normalizations obtained for the BFKL models are a factor of 6-9 lower than the 
data. 

Note that the BFKL 2^2 color-singlet fraction falls with A77 while the jet-level fraction 
rises. Although the BFKL cross section prediction has the same behavior with Ar] at both 
levels, parton showering for standard QCD processes reduces the jet A77 separation at the 
jet level relative to the parton level. The resulting BFKL color-singlet fraction prediction is 
thus sensitive to higher-order effects in QCD processes, casting doubt on the reliability of 
current parton-level predictions. 

Although the rate of t-channel photon exchange is too small to describe D0 data [[|, 
the photon-exchange process in HERWIG can be used to investigate two models in which the 
color-singlet couples only to quarks: a massless photon-like singlet and a massive U(l) gauge 



boson [24]. To simulate the exchange of a massive U(l) gauge boson that couples to baryon 
number, the photon cross section is multiplied by a kinematical factor of (1/(1 + m 2 B /E T )) 2 
which accounts for the exchange of a boson with mass tub- The coupling of the U(l) 
gauge boson is absorbed in the overall normalization. The behavior of the resulting color- 
singlet fraction is similar to photon exchange except for a steeper rise as a function of Et- 
The lowest allowed mass of m B = 20 GeV/c 2 [^4| gives the best fit to the data. Figure || 
shows the fits to the photon (dot-dashed line) and U(l) (dotted line) models. These models 
do not give an adequate description of the data due to their prediction of a steep rise in the 



color-singlet fraction with Et and with Ar] for the high Et sample [^ 

In addition to the color-singlet processes in HERWIG, we construct a set of simple "color- 
factor" models which explore different possibilities for the coupling of the color singlet to 
quarks and gluons. In this approach, the color-singlet fraction is parametrized as a weighted 
sum of the relative fractions of quarks and gluons from parton distribution functions. The 
weights represent effective "color factors" that enhance or suppress coupling to initial quarks 
and gluons compared to pure single-gluon exchange. This parametrization can be written 
as 



fs — In ■ (C qq ■ q P q P + Cqg ■ (q P g P + g P q P ) + C gg ■ g p gp) 

where /jy is a normalization factor, and qi and denote the relative fractions of initial quarks 
and gluons from hadron i (i — p,p). C qq , C qg , and C gg are the effective color factors for 
quark-quark, quark-gluon, and gluon-gluon processes, respectively |2B|, where the subscripts 
refer to the identity of the interacting partons from the proton and anti-proton. We apply 
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FIG. 5. Fits of Monte Carlo models to the color-sins; let fraction (a) ff D) {E T ) and (b)-(c) 
f^ D \Ari) for the low- Et sample and the high-Sy sample, respectively. Shown are comparisons to 
BFKL jet level (solid line), BFKL 2^2 parton level (dashed line), photon (dot-dashed line), and 
U(l) (dotted line) models. 



this formalism to several color-singlet models and make quantitative comparisons to the 
data. 

The soft-color rearrangement model, which has been used to explain rapidity gaps 
in positron-proton collisions ||15||, has recently been generalized as an alternative QCD- 



motivated explanation for rapidity gaps in hadron-hadron collisions fl6| . In this model, 
color flow (via the exchange of a single gluon) can be cancelled by the exchange of soft glu- 
ons, leading to a rapidity gap and thus an effective colorless exchange. This cancellation is 
more likely for initial-state quark processes than for gluon processes since quark states have 
fewer possible color flow configurations. Based on counting arguments with three quarks and 
eight gluons, a reasonable choice of color factors is C qq = l/9, C qg = l/2A and C 99 = l/64 |Jl6f . 
The "soft-color" model predicts a color-singlet fraction that increases with increasing par- 
ton x, and thus predicts a color-singlet fraction that increases with Et and Ar]. As seen in 
Fig. [| the soft-color model (dashed line) gives a reasonable description of the data. 

The exchange of two gluons in a color-singlet state was originally proposed as a simple 
mechanism to produce rapidity gaps between jets [0,0] with a predicted color-singlet frac- 
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tion on the order of 1% |TT|1. Within the naive color-factor formalism, 



for this simple two-gluon model are C qq = 1, C qg = 9/4, and C gg = (9/4) 2 |TT 



the color factors 
Since the 



coupling to gluons is stronger for the two-gluon singlet than for single-gluon exchange, the 
observed color-singlet fraction is expected to decrease with increasing x due to the parton 
distributions, in contrast with the soft-color model. As seen in Fig. ^|, the "two-gluon" model 
(dot-dashed line) does not give an adequate description of the data. 

We also fit a "single-gluon" model in which the color singlet is indistinguishable from 
single-gluon exchange, except in production rate. In this model, C qq = C qg = C gg = 1 and 
there is no dependence on or Arj (f$ is constant). The comparison of the single-gluon 
model and data is shown in Fig. || The single-gluon model (dotted line) does not reproduce 
the rising trend of the data with Et or Ar). 
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FIG. 6. Fits of color-factor models to the color-singlet fraction (a) J^ D \Et) and (b)-(c) 
f^ D \Arj) for the \ow-Et sample and the high-i?y sample, respectively. Shown are comparison 
to free-factor (solid line), soft-color (dashed line), single-gluon (dotted line), and simple two-gluon 
(dot-dashed line) models. 



To determine the best fit within the color-factor formalism, we fit a "free-factor" model 
in which the three color factors are allowed to float (/# = 1). We obtain the following values 
for the color factors: C qq = 0.025lg;ggJ, C qg = 0.00llg;gg?, and C gg = O^' 007 '. The best fit to 
the data, shown as the solid line in Fig. ^, thus favors a color-singlet that couples strongly 
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to quarks with only a few per cent contribution from the quark-gluon term (C qq : C qg : C gg = 
1:0.04:0). 

Table [D] shows the relative likelihood for each of the color-singlet models, ordered by 
relative likelihood. The relative likelihood is the likelihood of each model normalized to that 
of the free-factor model, which provides an excellent description of the data. Table [TT] also 
displays the average value of x 2 divided by the number of degrees of freedom ((x 2 )/df), which 
is calculated by integrating all possible \ 2 values over the prior probability distributions and 
dividing by the number of data points minus the number of fit parameters. If the Monte 
Carlo expectation were identical to the observed data in each bin then by definition (x 2 ) = 0. 

The soft-color model gives a good descriptions of both data variables, although the free- 
factor model has a likelihood that is several times larger than the soft-color model. The color 
factors for both of these models indicate that the data are consistent with a color-singlet 
that couples primarily to quarks. The single-gluon model does not give as good a fit to the 
data as the soft-color model, but a color-singlet that couples like a single gluon cannot be 
excluded. None of the other models adequately describes the color-singlet fraction for both 
variables. 



Model 


Fit toff D) {E T ) 


Fit to f {2D) (A V ) 


Likelihood 


(x 2 )/df 


Likelihood 


(x 2 )/df 


"free-factor" 


EE 1.0 


1/5 


EE 1.0 


6/9 


"soft-color" 


1/2.3 


3/7 


1/6.0 


10/11 


"single-gluon" 


1/28 


8/7 


1/350 


18/11 


BFKL 


1/10 5 


23/7 


1/66 


15/11 


photon 


1/660 


14/7 


1/10 5 


30/11 


"two-gluon" 


1/1100 


15/7 


1/10 5 


30/11 


BFKL 2^2 


1/10 8 


38/7 


1/10 9 


51/11 


U(l) 


1/10 14 


64/7 


1/10 11 


58/11 



TABLE II. Results of fitting the various color-singlet models to the color-singlet fraction 
measurements /^ D \Et) and fj? D \Arj). For each fit, the relative likelihood compared to the 
"free-factor" model and the (x 2 )/df (average \ 2 P er degree of freedom) are shown. 



Since the free-factor and soft-color models are in good agreement with the and Ai] 
dependence at y/s = 1800 GeV, we use them to calculate -Rigo . This results in a predicted 
value of -Rigoo = 1.5 ± 0.1, where the uncertainty includes the statistical error on the Monte 
Carlo samples as well as the difference between the two models. If we assume that the models 
are correct, then the discrepancy between the predicted value and the measured value of 
i?i800 = 3. 4 ±1.2 could be attributed to the a/s dependence of the survival probability. With 
this assumption, we obtain a value of 2.2 ±0.8 for the ratio of survival probabilities between 
the two energies, which is consistent with a recent estimate for this ratio of 2.2 ± 0.2 |18| . 

In conclusion, we have presented new information on the fraction of dijet events produced 
via color-singlet exchange. For our largest data sample, E T2 > 30 GeV at y/s = 1800 GeV, we 
measure a color-singlet fraction fs = [0.94 ± 0.04(stat)±0.12(syst)]%. The measured ratio 
of color-singlet fractions between i/i = 630 GeV and 1800 GeV for jets with Et2 > 12 GeV 
is 3.4 ± 1.2. The measured color-singlet fraction at = 1800 GeV tends to increase with 
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dijet Et and A77. Under the assumption of an ^-independent survival probability, the data 
favor soft-color rearrangement models in which the color-singlet is preferentially exchanged 
in processes with initial quarks. A model in which the color-singlet coupling is similar to 
single-gluon exchange cannot be excluded. Photon-like, U(l), and current two-gluon models 
do not adequately describe the data. 
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